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Abstract. We investigate whether the predictions of single-field models of inflation are ro¬ 
bust under the introduction of additional scalar degrees of freedom, and whether these extra 
fields change the potentials for which the data show the strongest preference. We study the 
situation where an extra light scalar field contributes both to the total curvature perturba¬ 
tions and to the reheating kinematic properties. Ten reheating scenarios are identified, and 
all necessary formulas allowing a systematic computation of the predictions for this class of 
models are derived. They are implemented in the public library ASPIC, which contains more 
than 75 single-field potentials. This paves the way for a forthcoming full Bayesian analysis 
of the problem. A few representative examples are displayed and discussed. 
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1 Introduction and Summary of the Main Results 

Inflation is the leading paradigm for explaining the physical conditions that prevailed in 
the very early Universe [1—6]. It consists of a phase of accelerated expansion that solves 
the standard hot big bang model problems, and provides a causal mechanism for generating 
inhonrogeneities on cosmological scales [7-12], The most recent Planck measurements [13-15] 
of the Cosmic Microwave Background (CMB) indicate that these cosmological perturbations 
are almost scale invariant, with no evidence of non-Gaussianities or isocurvature components. 
At this stage, the full set of observations can therefore be accounted for [16] in a minimal 
setup, where inflation is driven by a single scalar inflaton field 4> with canonical kinetic term, 
minimally coupled to gravity, and evolving in a flat potential V(<f)) in the slow-roll regime. 
When testing these models one by one in a systematic and Bayesian way [17-19], the best 
potentials, mostly of the “Plateau” type, can be identified. 

From the theoretical point of view, inflation proceeds at energy scales where particle 
physics is not known and has not been tested in accelerators. The physical details of how 
the inflaton is connected with the standard model of particle physics and its extensions 
therefore remain elusive. In particular, most physical setups that have been proposed to 
embed inflation contain extra scalar fields that can play a role either during inflation or 
afterwards. A natural question is therefore whether single-field model predictions are robust 
under the introduction of these extra fields. 

From the observational point of view, one should also focus on the extensions to the 
minimal scenarios that the coming generation of cosmological and astrophysical surveys will 
be able to rule out. As far as non-Gaussianities are concerned for example, the accuracy of 
currently proposed experiments is typically of the order [20-25] ct(/nl) ~ 1. Since models 
predicting |/nl| 1 are already excluded by Planck while models yielding |/nl| 1 will not 

be constrained in the short future, models that should draw our attention are those giving 

|/nl| ~ 1- 

For these reasons, in this paper we begin a systematic analysis of single-field slow-roll 
models of inflation when an extra light massive scalar field a is introduced and can play a role 
both during inflation and afterwards. 1 The potentials under scrutiny are thus of the form 
V{<j)) + m 2 cr 2 /2, where a is taken to be lighter than cj) at the end of inflation. Both fields are 
assumed to be slowly rolling during inflation, and eventually decay into radiation fluids with 
decay rates denoted T,* and T a respectively, during reheating. In the limit where the extra 
field a is entirely responsible for the observed primordial curvature perturbations, the class 
of models studied here is essentially the curvaton scenario [29-32], for which a natural value 
of the local non-Gaussianity parameter is [33] /nl = —5/4, hence of order one. 

However, in the present work, we address the generic setup where both cj) and a can a 
priori contribute to curvature perturbations [34-37]. When V{4>) is quadratic, these “mixed” 
scenarios have recently been studied in Refs. [38-40] and it has been shown [41] that the fit 
of quartic chaotic inflation can be significantly improved in the curvaton limit. In Ref. [42], 
a Bayesian analysis was carried out for the quadratic inflaton + curvaton models assuming 
instantaneous reheating, and these models were found not to be disfavored with respect to 
standard quadratic inflation. In this paper however, we aim at addressing all single-field 
potentials as recently mapped in Ref. [17]. Importantly, we do not make any assumption 

1 In many (notably string theory) high energy setups, a larger number of light scalar degrees of freedom 

are usually to be considered [26-28]. Here, for simplicity, we only implement a single additional light scalar 
field but the method we present can be extended to such scenarios. 
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as to the ordering of the three events: a becomes massive, cj) decays and a decays. Nor 
do we restrict the epochs during which a can dominate the energy content of the Universe. 
Assuming that the inflaton energy density decreases as matter after inflation and before it 
decays into radiation, this leaves us with 10 possible cases (including situations where a 
drives a secondary phase of inflation [35, 43-45]), see Fig. 2, that we study one by one. 
Furthermore, reheating kinematic effects are properly taken into account so that the number 
of e-folds elapsed between the Hubble exit time of the CMB pivot scale and the end of 
inflation is not a free parameter but is given by an explicit function of the potential U(</>) 
parameters, m a , r CT and cr en d, the vev of a at the end of inflation. This effect is particularly 
important for curvaton-like scenarios, since in these cases the same parameters determine the 
statistical properties of perturbations and the kinematics of reheating. It is therefore crucial 
to properly account for the interplay between these two physical effects and the suppression 
of degeneracies it yields. 

In practice, we make use of the SN formalism to relate observables to variations in the 
energy densities of both fields at decay time of the last field. This allows us to calculate 
all relevant physical quantities by only keeping track of the background energy densities. 
The main steps employed are described in section 2, while detailed formulas are provided 
in appendices A and B. The corresponding numerical routines have been implemented in 
the publicly available ASPIC code [46]. In this manner, one can now easily compute the 
predictions of the ~ 75 inflationary potentials originally contained in this library, when a 
massive scalar field is added, in all 10 cases mentioned above. 

In this paper, we discuss the results only for a few representative models: Large Field 
Inflation, Higgs Inflation (the Starobinsky model) and Natural Inflation. In appendix C we 
display their predictions in all 10 cases, in terms of scalar spectral index n s , tensor-to-scalar 
ratio r, and local non-Gaussianity parameter /nl- In general, isocurvature perturbations can 
also be produced [32, 47]. Their amplitudes depend on the details of the couplings between 
the fields 0 and a and the other species that here we leave unspecified. Indeed, they do 
not survive if either the field that last decays dominates the energy budget before decaying 
and/or the decay products fully thermalize after the last field decay. In principle however, 
non-adiabatic perturbations can give rise to extra constraints [48], and we leave this for future 
work. 

We describe and analyze our main results in section 3. However, in order to provide a 
quick picture of the situation, in Fig. 1, the predictions for n s and r are superimposed for 
the models mentioned above, when reheating is of the 7 th (left panel) and 9 th (right panel) 
types (see Fig. 2): 

• In the 7 th type, the inflaton decays before the a field becomes massive and then 
decays. This can be viewed as a standard “curvaton” scenario as originally proposed [29-32] . 
The tensor-to-scalar ratio r is driven to zero when the relative contribution from a to scalar 
perturbations increases. For small-field models (Natural Inflation, Higgs Inflation or Power 
Law Inflation in the figure), this does not make a big difference since those models already 
predict a value for r that is currently unconstrained. For large-field models however (Large 
Field Inflation in the figure), the observational upper bound on r strongly favors the model 
parameters for which r is reduced. On the other hand, the spectral index n s increases as 
the extra field contribution to the scalar power spectrum increases, reaching the asymptotic 
value n s ~ 1 — 2ei*, where ei* is the first slow-roll parameter evaluated at Hubble exit time 
of the pivot scale. For small field models, one effectively has n s ~ 1 in this limit and three 
cases can therefore be distinguished. If the single-field version of the model predicts a too 
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Figure 1. Predictions of a few inflationary models (Large Field Inflation V oc <j) p with p = 2, 4 
and 6, Higgs Inflation V oc [1 — exp(—-^/2/3^/Mp,)] 2 - the Starobinsky model Natural Inflation 
V oc 1 + cos(</>//) with / = 3 M P1 and Power Law Inflation V oc exp(—a<^/M P1 ) with a = 0.02), when 
a massive field a is added, and when the reheating scenario is of the 7 th (left panel, corresponding to 
the standard “curvaton scenario”) and 9 th (right panel, corresponding to a typical situation where the 
added massive field drives a secondary phase of inflation) type. In the left panel, the color encodes 
the relative contribution from a to the total curvature power spectrum. In the right panel, the 
color encodes the number of e-folds realized during the second phase of inflation, taking place during 
reheating. The black lines are the one and two sigma Planck 2015 contours [14], and the black squares 
stand for the predictions of the single-field versions of the models. 


red value for n s (for example, Natural Inflation in the figure), a small fraction of the mixed 
scenarios can provide a good fit to the data, otherwise n s is either too red or too blue. If the 
single-field version of the model already predicts the right value for n s (for example, Higgs 
Inflation in the figure), adding a massive field can only make the fit to the data worse. The 
same conclusion holds if the single-field version of the model predicts a too blue value for 
n s (for example, Power Law Inflation in the figure). For Large Field Inflation on the other 
hand, in the limit where only the extra field contributes to scalar perturbations, one has 
n s ~ 1 — 2 p/\p + 4(./Vend ~ -/V*)]. This means that p = 6 predicts values of n s that are too red 
and p = 2 predicts values of n s that are too blue in this limit, while p = 4 provides a good 
fit. 

• In the 9 th type (right panel of Fig. 1), a drives a secondary phase of inflation. For 
small-field models, the main effect is that lV enc i — IV*, the number of e-folds elapsed between 
the Hubble exit time of the pivot scale and the end of the first phase of inflation, decreases. 
(In particular, one can check that Power Law Inflation, which has an inflaton shift symmetry, 
has its predictions essentially unchanged). As a consequence, predictions are computed closer 
to the ending point of the first inflationary phase, where the potential V((j)) is less flat. This 
is why r typically increases and n s is shifted away from scale invariance. The situation is 
therefore opposite to the 7 th kind displayed in the left panel: models predicting too red 
values of n s provide even worse fit to the data once a is included, models predicting the 
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right value of n s are also made worse, while models yielding too blue values of n s (but still 
red) can provide a better fit for some parameters. For large-field models, this effect almost 
compensates with the reduction in r arising from the increase in the contribution from a to 
the total perturbations, so that r does not change much. In particular, r does not decrease 
sufficiently so that, even if one starts from a large single-field model with a too blue value of 
n s , a good fit to the data can never be obtained. 

In order to further investigate how the observational status of a given single-field model 
of inflation changes when a massive scalar field is added, one needs to incorporate non- 
Gaussianities to the discussion. This is done in detail in the rest of the paper. The aim of 
the above discussion is only to provide the reader with a taste of the physical effects we will 
encounter later. Finally, let us note that in some reheating scenarios and for some potentials, 
the best fit to the data can be improved when adding a massive scalar field, but this might 
be realized only for limited ranges of parameters. The immediate question to ask is of course 
whether the increase in the fit and the corresponding volume in parameter space could lead to 
a substantial increase in the Bayesian evidence of such models. In a future publication [49], 
we will address this issue by computing the Bayesian evidence of the models present in the 
ASPIC library [46], when a massive scalar field is added. This will allow us to quantitatively 
discuss the observational status of curvaton-type models of inflation, as well as assessing the 
robustness of single-field models of inflation when massive scalar fields are added. 

2 Method 

In this paper, we investigate the situation where inflation is mainly driven by an inflaton scalar 
field (f) with potential V(</>), in the presence of a light free scalar field a. Here, “light” must 
be understood in the sense that the mass of this extra field is smaller than the Hubble scale 
during the whole inflationary period, i.e. m a /V 2 <C H c n d, where H e n d denotes the Hubble 
factor at the end of inflation. The two fields follow Klein-Gordon equations in Friedmann- 
Lemaitre-Robertson-Walker space-times given by 

(f> + 3 H<j) + = 0 , 

a + 3 Ha + m 2 cr = 0 , 

where the Hubble factor H obeys the Friedmann equation 

3M P 2 ,F 2 = y + V (4) + y + ^fa 2 . (2.3) 

During inflation, we assume that both fields are slowly rolling. When inflation is over, 
we work under the sudden decay approximation, since it has been shown to be a good 
approximation to the full numerical results in the usual curvaton scenarios [50, 51]. In 
practice, this means that when H drops below (respectively T^), the energy contained 
in cf> (respectively a) instantaneously decays to radiation. 2 We also assume that close to its 
minimum, the inflationary potential is quadratic so that the energy density contained in <j) 
behaves as matter between the end of inflation and the decay of cj>. 

2 Here, F^, (respectively r CT ) are effective values for which assuming instantaneous decay at H = Y ^ (re¬ 
spectively H = T a ) provides a good description of the full decay dynamics. They do not exactly match the 
“decay rates” as would be defined in a dynamical way. However, since T $ (respectively r CT ) are varied across 
orders of magnitudes, this plays a minor role. 


( 2 . 1 ) 

( 2 . 2 ) 
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Figure 2. Different possible reheating scenarios, depending on the values taken by T CT , rn a , T^, H en d 
and a en d- Cases 1, 2 and 3 correspond to T^, < r CT < m a < H e n( j; cases 4, 5 and 6 correspond to 
r CT < < m a < H e n( j; cases 7, 8, 9 and 10 correspond to r CT < m a < < H en d. Within each 
row, different cases are distinguished by a en d/M P1 which controls when a dominates the total energy 
density. The blue curves stand for the energy density of <p while the green ones are for a. 


On top of the parameters describing the inflationary potential V((/)), the system therefore 
contains four free parameters: m a , and a em j. The kinematic parameters of reheating 

(energy density at the onset of the radiation era and mean equation of state parameter, see 
section 2.3) being entirely fixed by these quantities, there is no other free parameter. 

Let us now see what the different scenarios for the reheating phase are. Here, by 
“reheating”, we refer to anything happening between the end of the phase of inflation driven 
by 4> and the onset of the final radiation dominated era. Four energy scales are present in 
the model, r CT , m a , and H e n d- Since m a < H en( i by definition, and since m a > T a (i.e. 
we assume perturbative decay) and H e n( j > (0 decays after inflation ends), these values 

must be ordered according to one of the three following possibilities: < m a < H e n d, 
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r fj < r 0 < m a < H cn d or To- < m a < T^ < H e nc j. The time evolution of the energy densities 
associated with </> and a are schematically represented in Fig. 2. Depending on the value 
of cr en( j, hence on the energy density contained in a. one can see that 10 sub-cases must be 
distinguished. 


2.1 Observable Quantities as Energy Density Variations 

Since we want to address a large number of inflationary potentials V(cj)) and describe all 10 
possible cases in a systematic way, it is necessary to design a calculational framework that 
allows us to obtain the predictions of these models in a standardized and methodical way. A 
numerically-oriented well-optimized solution is to make use of the “quasi-isotropic” [52-54] 
or “separate universe” [55-58] approach and the corresponding 5N formalism [56, 59-61] to 
express all observable quantities in terms of the variations of the energy densities at the last 
decay time, with respect to the fields values at first Hubble crossing time during inflation. 
Let us see how this can be done. 

When the last field decays, since the other field has already decayed, one has two fluids: 
a matter fluid that corresponds to the field ip m oscillating at the bottom of its potential, and 
a radiation fluid that corresponds to the decay products of the other held <p r . In cases 1-3, 
one has <p m = 4> and <p r = cr since the held that decays last is </>, while in cases 4-10, one has 
<p m = a and <p r = (j) since the held that last decays is a. 

The non-linear curvature perturbation of the matter fluid on uniform-matter density 
hypersurfaces is given by [58] 


Cm {t, x) 


SN (t,x) + f 
J p 


Pm(t) 


d/0 m 
3 Pm 


C {t, x) + - In 


Pm(t,x) ~ 

Pm(t) J 


(2.4) 


where in the last equality p m (t,x) is defined on uniform total energy density hypersurfaces. 
In the same manner, the non-linear curvature perturbation of the radiation fluid on uniform- 
radiation density hypersurfaces is given by, 


Cr (t, x) 


C {t, x) + - In 


' Pr(t,xY 
. pt(t) \ 


(2.5) 


where similarly, p T (t,x) is dehned on uniform total energy density hypersurfaces. The total 
density is uniform on the last decay surface, parametrized by H = T, so one has p m (tdeo x ) + 
Pr (f-dec i x) = p(idec)> where a bar denotes background, homogeneous quantities. This gives 
rise to [51] 

Pm (idee) e 3(<:m_C) + p r (t dec ) e 4(Cr-C) = Pm (idee) + Pr (W) • (2-6) 

This non-linear relation allows us to relate C with Cm an d Cr at an y order in perturbation 
theory. 


2.1.1 Power Spectra 

Let us now expand Eq. (2.6) at linear order to work out the scalar power spectrum. Defining 


Liec 


3 P m (td e 


3p m (td ec ) T 4p r ( t d ec) 


(2.7) 


one obtains C = r decCm + (1 — r^ec) Cr- We thus need to specify Cm and Cr- Since the two 
corresponding calculations proceed in exactly the same way, let us give the details only for 
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Cm- A first remark is that, from Eq. (2.4), one can write 


/- _ 1 ! Pva.it , x) 

C " “ 3 ln [ *„(«) 


8N=0 


( 2 . 8 ) 


On the other hand, after Fourier expanding the density field p m , one has, on super- 
Hubble scales, 


Spm(t, |,5JV=0 


dp m (t) dp m {t) 

TT7- Op m ,fc + -7T7- 

5tv=0 


8N=0 


(2.9) 


Here, Sp mt k and <5(/j rj fc are the held fluctuations at the Hubble exit time of the wavenumber 
k. At leading order in slow roll, they are uncorrelated Gaussian fluctuations with amplitude 
H{k)/{ 2ir), where H(k) is the value of the Hubble parameter when k exits the Hubble radius 
(see Ref. [62] for inclusion of higher order slow-roll corrections). Here also, |(5iv=o means that 
the derivatives with respect to pk, the values of the fields when k exits the Hubble radius, 
are to be evaluated under the constraint that the number of e-folds elapsed between the time 
IVfc when the mode k crosses the Hubble radius during inflation and the time N^ ec of the last 
decay, is fixed. Plugging Eq. (2.9) into Eq. (2.8), one obtains, at linear order, 


Cm (t,k) = 


3pm(t) _ d(/) m ,k 


dp m (t) 


(*) 

OPm,k + o / 
SN=0 94>T,k 


fi(Pv,k 


8N=0 


( 2 . 10 ) 


A similar expression for Cr can be derived, where one just needs to exchange matter and 
radiation indices, and 1/4 replaces the overall 1/3 factor. One then obtains, at the last field 
decay time, 


C (k) = 


^dec 

dPm(tdec) 

+ 

fi(N dec —N k )=0 

1 ^dec 

dprit dec) 


3/Jm(^dec) 


4pr(^dec) 

9(f>m ,k 

S(N iec —N k )=0_ 

^dec 

dpmitdec) 

+ 

fi( N dec— N k)=0 

1 - ^dec 

^Pr(^dec) 


3pm(^dec) 

,k 

4pr(^dec) 

,k 

$(N dec —N k )= 0_ 


+ 


In order to make notations more compact, let us introduce the quantities 


A% (A;) = ain f a (tdec) 
p dpp,k 


$Pm ,k 

$Pt ,k • 

( 2 . 11 ) 

( 2 . 12 ) 


<5(A r dec-A r i;) = 0 


where a,/3 = nr,r. The scalar curvature power spectrum is then given by V^(k) = (C 2 (&))j 
that is to say 


V C (k) = 


r dec i 1 r dec , r 


- A m 4- 

3 Am + 


4 Am I + 


r dec A m , 1 r dec , r 


-A™ + 
3 r 


4 * 


HJM) 

2vr 


-i 2 


(2.13) 


From here, the spectral index and the running can be obtained as simple expressions of the 
potential derivatives at first [63] and second [62] order in slow roll. For simplicity, here we do 
not reproduce these formulas since they are explicit. From Eq. (2.13), one can also identify 
the relative contribution from the matter field to the total scalar power spectrum (first term 
in the brackets) or, equivalently, the relative contribution from the radiation field (second 
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term). These quantities are not directly observable by themselves but they can be useful to 
study as they indicate which field mostly contributes to scalar perturbations. 

Let us now consider the gravitational waves power spectrum. The tensor perturbations, 
in contrast to the scalar perturbations, remain frozen on large scales and decouple from 
the scalar perturbations at linear order. Thus the primordial perturbation spectrum for 
gravitational waves is given by the usual formula [64], 

-w-Kff 


2.1.2 Non-Gaussianities 

Due to the presence of non-adiabatic field perturbations on large scales after inflation, 
curvaton-like scenarios can produce local type non-Gaussianities [51], which is also the ob- 
servationally most constrained type of non-Gaussianities [15]. This is why in this paper, we 
characterize non-Gaussianities predictions by the value of the local /nl parameter. It can 
be calculated by expanding perturbations up to second order, C = Cl + C 2 / 2 , where Cl is the 
quantity we already calculated for the scalar power spectrum in Eq. (2.11) in section 2.1.1. 
The local non-Gaussianities parameter /nl is defined through the ratio of the bispectrum to 
the power spectrum squared, 


5 (C 3 ) 5 (Cfo) 

6 3(C 2 ) 2 6( C 2) 2 


(2.15) 


where the second expression is valid at leading order in the perturbations. Let us see how it 
can be calculated. 

Similarly to what was done for the scalar power spectrum, one can first expand Eq. (2.6) 
at second order in the perturbations to relate C 2 at the last decay time to its matter and 
radiation components. One obtains 


C 2 — LiecCm2 T (1 ^dec) Cr2 T l"dec (1 Llec) (3 T ^dec) (Cml Crl) • (2.16) 


The first order curvature perturbations Cml and Crl have already been expressed in terms of 
the scalar field fluctuations at Hubble exit time in section 2.1.1, see Eq. (2.10). One needs 
to work out similar expressions for Cm 2 and Cr 2 - Expanding Eqs. (2.4) and (2.5) at second 
order in the perturbations, one obtains 


1 d 2 In p m 

Cm2 («) = r 
6 


d(j) 2 


m ,k 


2 1 d 2 lnp m 

+ e ~a^~ 

SN=0 T ’ k 


SN=0 


and a similar expression for Cm 2 - Recalling that, 
K, fc > = K*> = H (k) 7(2vr) 2 , <*<*) = K fe > 



d 2 In p m 


^'r'n].k r.A- , 

SN=0 


(2T7) 

at Hubble crossing time, one has [65] 
= 2H(k)^/ (27t) 4 and (6<p m)k 5ip I}k ) = 0, 
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[1 + 6 (r dec - 1) r dec ] + rdec {27 (Ar) 2 (r^ + r dec - 5 ) - 24A™ m - 12A” 

+8 (A™) 2 [(6r dec (2 + r dec ) - 17] - 36A r m AJ [r dec (3 + 2r dec ) - 8]})] } / 

{9 [KJ 2 + K) 2 ] (r dec - l) 2 - 24 (r dec - 1) r dec (A£A^ + A r m A^) 


+16 


(Am ) 2 + (A™) 2 


r 2 

dec 


} 


(2.18) 


Here, similarly to Eq. (2.12), for display convenience, we have introduced the short notation 


A%{k) 


d In p a (i d ec) 


( 9+/3,/c^+7,fc 


H N dec—Nk)=0 


(2.19) 


This formalism might seem a little heavy (in particular, in light of Eq. (2.18) or consid¬ 
ering the various equations given in appendix B), but it allows us to design a computational 
strategy that is fully generic, since it only relies on tracking down energy density components 
during the different stages of the scenarios under consideration. When evaluated at the last 
field decay time indeed, these give rise to rd ec through Eq. (2.7). Then, fixing and N dec to 
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their background values, we just need to take the derivative of these quantities with respect 
to (j)k and <7fc up to second order in order to calculate A ^ and defined in Eqs. (2.12) 
and (2.19). Observable quantities are finally expressed in terms of these parameters only 
through Eqs. (2.13), (2.14) and (2.18). This is why in the next sections, we calculate how 
energy densities evolve throughout inflation and reheating. 


2.2 Inflation 


During inflation, we assume that both fields 4> and o are slowly rolling. In terms of the 
number of e-folds N = ln(a), Eqs. (2.2) and (2.3) then give rise to 


= _ YL _, 

dN P1 V + m^ff 2 /2 


do - , r n m 2 a 

dlV lvlpi V + m 2 a 2 /2 ’ 


( 2 . 20 ) 


where a prime denotes derivative with respect to 0. Combined together, these equations 
yield der/d0 = m 2 a/V'. This relation can be integrated, and one obtains 


O’ * 


(Tend exp 



d <f> 
V 7 


( 2 . 21 ) 


On the other hand, combining Eqs. (2.20) and (2.20) also gives rise to V/V'dcj) + cr/2da = 
—MpjdlV, so that one has 



^end — A* . 


( 2 . 22 ) 


Plugging Eq. (2.21) into this relation leads to the number of e-folds realized between and 

0end j 


iV en d - A* 





(2.23) 


When A^ en d — IV* is known, can then be computed (often numerically) from inverting 
Eq. (2.23) and cr* can then be derived from Eq. (2.21), provided </> en d is specified. This can 
be obtained from requiring that the first slow-roll parameter, 


ei ~ 


Ml 


V' 2 + mta 2 




2 ’ 


(2.24) 


equals one at the end of inflation. One then has <fi 2 /2 + a 2 / 2 = V((j))/2 + m 2 / 2, so that 
Pend = (3/2)E(0end) + 3m^(7g nd /4 an d ^cnd can be computed through the Friedmann equa¬ 
tion (2.3). From here, variations with respect to (ft* and <7* of cr en d, A en d and the energy 
densities contained in both fields at the end of inflation can be calculated, and we display 
the corresponding formulas in appendix A. 


2.3 Reheating 

During reheating, the way energy densities evolve is different for the ten cases under consid¬ 
eration. In appendix B, we detail the calculation of rd ec and the energy density variations 
case by case, and we derive under which conditions on the parameters m CT , r^, T^ and cr en d 
each scenario takes place. 
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It is also important to mention that the parameters of the model completely fix N end — 
N *, the number of e-folds elapsed between the Hubble exit time of the pivot scale k * around 
which observable quantities are expressed, and the end of inflation. This can be seen intro¬ 
ducing the reheating kinematic parameter, 


In R r a( j 


1 3u> re h j / Preh 

12 (1 T" Wreh) \Pend 


Here, the averaged equation of state w re h is defined as 


^reh 


1 

AiV re h 



w (A) d N . 


(2.25) 


(2.26) 


As explained in Refs. [66, 67], the reheating parameter is related to N m d — IV* through 


A en d - A* 


In R iad T In 


(1 ft 

V 27 MpiPend 1 


1 

ei*/3 


— In 


k * / ttn 


-1/4 
Pi'.now 


(2.27) 


where p r ,now is the energy density of radiation today rescaled by the number of relativistic 
degrees of freedom. The reason why In i ? rad is convenient is that it follows a simple additivity 
rule. 3 In particular, this means that any radiation dominated phase will not contribute to 
lnR rac j (since, in this case, 1 — 3tc re h = 0). This is consistent with the fact that reheating is, 
by definition, the connection between the end of inflation and the beginning of the radiation 
era. 

In all ten cases detailed in appendix B, reheating is comprised of three kinds of sub¬ 
phases: matter phase, radiation phase, and inflating phase. For matter dominated phases 
occurring between p\ and p 2 , w = 0 and one simply has 


ln Kid = — In 


12 



(2.29) 


When an extra phase of inflation occurs during AIV in f e-folds and proceeds between p\ and 
P 2 , one has 4 

IniCi = -AAT inf + 1 ln . (2.31) 

Finally, as already noticed, during a radiation phase, ln 7?, t r . = 0. For each reheating scenario, 
one then just add to sum the contributions of all sub-phases to compute lni? ra d, hence 

A en d — A*. 

3 Indeed, let us consider the case where reheating is made of two phases. The first one drives the energy 
density from po to pi with an average equation of state rhi, and lasts iVi = ln(po/pi)/[3(l + wi)] e-folds. 

Similarly, the second phase drives the energy density from pi to p 2 with an average equation of state W 2 , 
and lasts N 2 = ln(pi/p2)/[3(l + 102 )] e-folds. From here, one can calculate the averaged equation of state 
parameter iZ} re h = (JViwl + N 2 W 2 )/(Ni + N 2 ), and the reheating parameter 


ln R ra a — 


1 — 3wi 
12{l + m) 


hr 


+ 


1 — 3u’2 
12(1 +w 2 ) 


ln 


1 3rC r eh 
12 (1 + w leh ) 


which is simply given by the sum of the reheating parameters during each phase. 
4 The integral in Eq. (2.26) can be expressed as 


/”“ = /(- 


d H 


H 


ln 


-1 + ^ei ) AN = -AN- Z I = -AN + i hi ( — 


P 2 


(2.28) 


(2.30) 


where we have used the Friedmann equation in the last equality. This gives rise tow = — 1+ln (pi/p 2 ) / (3AN). 
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3 Classification and Discussion 


The numerical routines corresponding to the procedure detailed in section 2, with relevant 
formulas given in appendices A and C, have been added to the publicly available ASPIC 
library [46] and can be downloaded from the ASPIC website. This means that one can now 
easily compute the predictions of the ~ 75 inflationary potentials originally contained in this 
library, when a massive scalar field is added, in all 10 reheating scenarios. For illustrative 
purposes, in appendix C, we display these predictions for three prototypical examples: 

• Large-field inflation, for which V = AI 4 ((f>/M P1 ) p , is presented in appendix C.l. It is a 
typical example yielding a value for r which is too large in its single-field version. 

• Higgs inflation (the Starobinsky model) is presented in appendix C.2, and its potential 
is given by V = M 4 [ 1 — exp(—y / 2/30/M P1 )] 2 . It is a typical example for which the 
single-field version of the model provides a good fit to the data, both in terms of n s 
and r. 

• Natural inflation is presented in appendix C.3, and its potential is given by V = M 4 [l + 
cos(</>//)]. When / is not super-Planckian, it is a typical example that yields a value 
of n s which is too small in the single-field version of the model. 

In all 10 cases, predictions are plotted in the (n s ,r), (ra s ,/NL,) and (/nlU’) planes. They 
are compared to the predictions of the single-field versions of the models, if one takes the 
averaged equation of state parameter for the reheating epoch of these models to be w ie h = 0. 
Indeed, in the present work, let us recall that we assume that close to its minimum, U(</>) 
is quadratic so that the energy density contained in <j) behaves as matter between the end 
of inflation and the decay of 4>. For fair comparison purpose, we therefore make the same 
assumption regarding the single-field versions of the models. Finally, the Planck observational 
constraints are also superimposed. 

Let us describe the behaviors obtained in all ten cases. In order to gain some intuition 
on the order of magnitude of the predicted amount of non-Gaussianities, a useful formula is 
the estimate [41, 47, 51] 

/NL ~il(n) (i" 4_2rt, “)’ (3,1) 

where V?/Vq is the fractional contribution from a to the total curvature power spectrum, 
defined below Eq. (2.13). In particular, one can see that the only way to get large non- 
Gaussianities, is to have a providing a substantial contribution to the total amount of scalar 
perturbations while the field that first decays remains dominant afterwards (so that r^ ec <C 1). 
We now review the different cases, and schematically summarize the main observed trends 
in table 1. 

Case 1 is in practice indistinguishable from the single-field version of the model since 
c7 never dominates the energy content of the Universe and never provides the dominant 
contribution to the scalar perturbations. It decays before 4> and is a purely spectator field. 

Case 2 is the only situation where a dominates the energy content of the Universe for 
a transient phase only. It even drives a short second phase of inflation. It decays before </> 
and is subdominant at the last field decay, which means that, again, it does not contribute 
significantly to scalar perturbations. This is why the predicted values for n s and r are 
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Reheating 

Extra 

r dec 

V C /'Pc^end) 

re s( cr end) 

r{cr end ) 

/nl 

Scenario 

Inflation 

SF 

LF 

SF 

LF 

SF 

LF 

1 

no 


o/o 

★ 

★ 

<i 

2 

yes 

~1 

o/o 

★ 

★ 

<0(0.1) 

3 

yes 

<1 

0^0 

0 1 

* 

\ 

* 

-k 

<1 
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no 

<1 

0 /* 1 

* 1 

* 1—2ei» 

* s \i 0 

>1 

5 

no 

~1 

1 \ 0 

1 * 

1—2ei* * 

0 Z * 

0(1) 

6 

yes 

~1 

0 /• 0 

0 1 

* 

\ 

* 

■k 

<1 

7 

no 

<1 

0 1 


* Z 1—2ei* 

* \ 0 

>1 

8 

no 

~1 

1 \ 0 

1 N \/ ★ 

1—2ei* * 

o 

\ 

* 

0(1) 

9 

yes 

~1 

0^0 

0 1 

* 

\ 

* 

■k 

<1 

10 

yes 

~1 

0 Z' 0 

0 1 

* 

\ 

* 

■k 

-Cl 


*: single-field predictions 


Table 1 . Schematic trends observed for the 10 reheating scenarios. The second column shows whether 
an extra phase of inflation takes place during reheating or not. The schematic value given for the 
ratio rdec defined in Eq. (2.7) indicates whether the matter field dominates the energy content of the 
Universe at the last field decay time (rd ec ~ 1) or whether the radiation field does (rd ec < 1). The 
ratio T° /7\- shows how the relative contribution from a to the total scalar power spectrum evolves 
with (Tend) and between which typical values it varies. In some cases, the trends depends on whether 
the single-field version of the model under consideration is of the small-field (“SF”) or large-field 
(“LF”) type (where small and large field must be understood in terms of ei*). The same trends are 
shown for n s and r. Here, the symbol * refers to the single-field version of the model predictions. 
Finally, we report whether large non-Gaussianities can be produced. 


again almost the same as the ones for the single-field version of the model. However, non- 
Gaussianities are slightly enhanced, and can typically reach up to 0(0.1) values (which, 
though small, are significantly larger than the single-field predictions). 

Cases 3, 6, 9 and 10 have similar phenomenologies. In these cases, a secondary phase of 
inflation takes place during reheating. The difference with case 2 is that this phase does not 
necessarily need to be short, and importantly, at the last field decay time, a dominates the 
energy budget of the Universe. This explains why non-gaussianities remain small. However, 
the behavior of n s and r depends on whether the single-field version of the model is of the 
small or large field type. For small-field models, the main effect of increasing <7 en d is that 
more e-folds are realized in the second phase of inflation. As a consequence, A r enc ] — N * 
decreases, and the predictions are calculated at a field location <^>* that gets closer and closer 
to <p e nd, where the potential is steeper. This means that n s is shifted away from 1 and r 
increases. Let us stress that the interplay between u en d, the number of e-folds realized in the 
second phase of inflation and IV enc j — IV*, is described by the reheating kinematic analysis of 
section 2.3. A proper inclusion of these effects is therefore crucial to correctly account for 
these regimes, where reheating kinematic effects dominate. When <7 en d continues to increase, 
the contribution from a to the total amount of scalar perturbations starts to increase as well. 
However, before it reaches non vanishingly small values, the number of e-folds realized in the 
second phase of inflation is too large so that the scales that we observe in the CMB would 
not have exited the Hubble radius during the first phase of inflation. For large-field versions 
of the single-field model, on the contrary, both effects have the same order of magnitude: 
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adding more e-folds during reheating drives </>* closer to (f> en d (hence n s is shifted away from 
one, as in the small-field case), but increasing (7 en d also enlarges the contribution from a to 
the total amount of scalar perturbations. As far as r is concerned, both effects act in opposite 
directions and roughly cancel out (especially if the potential is quadratic, where, up to mild 
logarithmic terms, they exactly cancel out), which implies that r is almost unchanged. 

Cases 4 and 7 give comparable phenomenologies. The extra field a remains subdomi¬ 
nant throughout the entire scenario but contrary to case 1, it decays after </, which means 
that as (Tend increases, its contribution to the total amount of scalar perturbations becomes 
significant. In this regime, large levels of non-Gaussianities can therefore be produced. As 
far as n s and r are concerned, one can show that they are well approximated by [63] 

n S ~ 1 ~ - p-J ("S - 1)single—field “ 2 D* , (3.2) 

f v c\ 

r ~ ^single—field ■ (3-3) 

When cr en d is small, represents only a small fraction of Vq and one can check that the 
same predictions as in the single-field version of the model are obtained. In the opposite 
regime, V° ~ r takes vanishingly small values and n s reaches 1 — 2ei*. For small-field 
models, this is essentially not distinguishable from n s = 1. For large-field models however, 
n s can be substantially different from 1 and one can see that, for example, V oc provides 
a good fit to the measured value of n s in this regime. 

Cases 5 and 8 are similar to these two previous cases except that a dominates the energy 
budget of the Universe when it decays. This implies that large non-Gaussianities cannot be 
produced since rd ec ~ 1- In fact, when <r en d is not too large, a provides the main contribution 
to the scalar perturbations and one has /nl ~ —5/4, see Eq. (3.1). When cr en d increases, its 
contribution to Vq decreases, /nl takes 0(1) values and eventually becomes much smaller 
than one. As far as n s and r are concerned, they interpolate between the same values as 
for cases 4 and 7. The only difference is that the relative contribution from a to V £ now 
decreases with cr en d, hence <r en d varies in the opposite direction between the two asymptotic 
values. 

4 Conclusion 

The ongoing quest for high accuracy astrophysical and cosmological data offers an unprece¬ 
dented opportunity to constrain the inflationary theory. So far, observations are compatible 
with the simplest framework, where inflation is driven by a single scalar held, minimally 
coupled to gravity and slowly rolling down a hat potential with a canonical kinetic term. 
However, most physical setups that have been proposed to embed inflation contain extra 
scalar degrees of freedom that can play a role either during inflation or afterwards. This is 
why it is important to ask the two following questions: Are single-held models predictions 
robust under the introduction of these extra helds? Do the data show preference for the same 
single-held potentials once these extra helds are included? 

This paper is the hrst of a series which tackles these two issues. It investigates the 
situation where an extra light scalar held is added to single-held scenarios, and contributes 
both to the total curvature perturbations and to the reheating kinematic properties. Since 
it has been shown [68] that the later plays a crucial role when constraining models from 
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measuring the former, it is indeed necessary to properly describe these two effects simul¬ 
taneously. In practice, ten different types of reheating scenarios have been identified and 
studied separately for each single-field potential. Making use of the 5N formalism to relate 
observable quantities to variations in the energy densities contained in both fields, we have 
derived all necessary formulas to implement a systematic and “industrial” treatment of this 
class of models. The corresponding routines have been added to the publicly available run¬ 
time library ASPIC [46]. It now allows the fast computation of the predictions of the ~ 75 
inflationary potentials originally contained, when an extra massive scalar held is added and 
in all 10 reheating scenarios. 

For illustrative purpose, we have displayed the predictions obtained for a few potentials. 
The few generic trends we have thus identified have a number of observational consequences 
for the models we considered: 

• In their single-field versions, large-field models are plagued with values of r that are 
too large. We saw that scenarios producing a second phase of inflation during reheating 
do not help to solve this issue. The only scenarios where r is decreased are 4, 5, 7 and 8. 
However, in cases 4 and 7, large non-Gaussianities are produced and one concludes that these 
models are safe only in scenarios 5 and 8. The question then is how much parameter space 
allows one to obtain the correct predictions and how much these models are improved. 

• Models that already provide a good fit to the data in their single-field version, such 
as Higgs Inflation (the Starobinsky model), can only perform worse once an extra scalar 
field is introduced and shifts the predictions away from the preferred values. The question is 
again how much parameter space is wasted in these disfavored regions and how it affects the 
likelihood of these models. 

• For small-field models producing values of n s that are too small, such as Natural 
Inflation, the only way to increase n s is to consider cases 4, 5, 7 or 8. However, cases 4 and 
7 yield large non-Gaussianities and one is left with cases 5 and 8. There, in the asymptotic 
regime, n s is driven to 1 so there exists only an intermediate range of values for the model 
parameters such that n s acquires the right amplitude. The question is of course how fine 
tuned this range of values is. 

To address these issues, one needs to resort to a complete Bayesian analysis of the 
problem. More generally, this is the only way to make accurate statements about the obser¬ 
vational status of the models we considered, the robustness of single-field models of inflation 
under the introduction of extra scalar degrees of freedom, and the constraining power of the 
current and future data to prove or disprove such models. Such an analysis will be presented 
in a separate publication [49]. 
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A Energy Density Variations During Inflation 


The slow-roll trajectory (2.21) gives rise to 
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(A.5) 


where the derivatives of AndAend) are given by varying </> end and cr en d in the condition 
eiend = 1 an d making use of Eq. (2.24). Combined with Eqs. (A.1)-(A.5), they also give rise 
to the derivatives of And (A V*)- Making use of the Friedmann equation H 2 nd = [2V ((/> en d) + 
3m 2 cr 2 nd ]/(4M 2 1 ), the derivatives of H cn d with respect to A and a* can also be obtained. 
Then, noting that 
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(A.6) 


the derivatives of p^ nd with respect to A and a* can be derived. We do not reproduce them 
here since they are straightforward. Finally, the number of inflationary e-folds (2.22) obeys 
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B Energy Density Variations During Reheating 

B.l r ^ <c r CT <c 777.(7 <C H end 

In this case, the curvaton becomes massive and decays before the inflaton decays. This 
corresponds to cases 1, 2 and 3. Since the last field to decay is the inflaton, one has tp m = 0 
and tp r = a. 


B.1.1 Case 1 


In this case, between the end of inflation and the decay of the inflaton, the Universe is 
dominated by the inflaton which behaves as matter, and one has 

H= - 3 Hend -. (B.l) 
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The curvaton field a becomes massive when H a _ mass = m a /\/2. Together with Eq. (B.l), 
this gives rise to 
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Assuming that a is slowly rolling between the end of inflation and the moment it becomes 
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Combined with Eq. (B.2), this gives rise to 
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One then has p °_ mas s = 3m 2 <7 2 _ mass /4. Then o simply behaves as matter, and one has 


Pa —dec Pa —r 


^a —dec 


CLa — 


a —mass 


-3 


°-n2 2 

= ^^end eX P 


4 

_ 9 + 9i7 2 


Per—mass 
2m 2 




dec 


ik- 


cr—mass 


end 


(B.5) 

(B.6) 


Then a decays and the associated energy density evolves as radiation. At time of inflaton 
decaying, one then has 
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As for the inflaton, since it dominates the energy content of the Universe when it decays, its 
energy density is simply given by 


P ™-dec — — p^_ dec , (B.9) 

and r dec can be computed making use of Eqs. (2.7), (B.8) and (B.9). 

Let us now work out the perturbed energy densities. They must be evaluated at the 
fixed background time A^_ dec , so remember that this quantity must remain unperturbed. 
For the matter component, one has 
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Let us now calculate the reheating parameter. Since the curvaton never dominates the 
energy budget of the Universe, and reheating is only made of a matter dominated phase 
taking place between H e nc j and T^, one simply has 


In R rad = - In 
o 



(B.16) 


Finally, let us determine under which condition case 1 takes place. In the limit where 
m a < H e nd , one has /?£_ mass ~ 3m^^ nd exp(-4/9)/4, while p^_ mass ^ 3M^ml. Requiring 
that p°_ mass < Pt -mass then leads to 


°end 

M P , 


« Vie 2 / 9 


(B.17) 


B.1.2 Case 2 

In this case, the curvaton comes to dominate the energy content of Universe while it is still 
slow rolling, and let the inflaton dominate again after a decays. Let us denote t eqi and f e q 2 
the two times when (f = p a . At time of first equality between p a and p one has 


p* = p* ( 3 = p<t> ( 2 

Pe qi P e nd{ aend ) ^V^cndJ ' 


(B.18) 


On the other hand, one can combine Eqs. (B.l) and (B.3) to get 
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which gives rise to 
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(B.20) 

The Hubble parameter H e q can then be obtained numerically equaling Eqs. (B.18) and (B.20), 
and er eqi follows from Eq. (B.19). Then follows a second phase of inflation driven by 
a, which ends when the curvaton becomes massive, that is when cr CT _ mass = \/2M P1 and 
H a -mass = m a /y/2. The number of e-folds realized in this phase is given by iV CT _ mass — N eqi = 
(<4n - mass)/( 4M li) = /( 4M ei) “ V 2 - As a consequence, we have 
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Pa —r 


= Pe qi exp [-3 (iV CT _ mass - N eqi )] = pi nd exp ^ 


3tJ e qi ' 

4Mp! 


(B.21) 


On the other hand, one simply has p£_ mass = 3/4m^cr^_ mass = 3/2rn^Mp V Then, the 
Universe is dominated by a matter fluid, until the curvaton decays. One then has 
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(B.22) 
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while the energy density of the curvaton is simply p a _ dec — SM^Tl. Then follows a phase 
where the Universe is dominated by the decay products of the curvaton, which behave as 
radiation, until the inflaton dominates the energy budget again. At this time, the energy 
density of the inflaton is given by 
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(B.23) 


while the energy density contained in radiation is p r cq2 — — ptq 2 - By equating these 

two quantities, one obtains 


H eq2 — 





(B.24) 


Then, the Universe is dominated by the inflaton again, which behaves as matter, until it 
decays. The energy densities of the matter and radiation fluids at this time are respectively 
given by 
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and 
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2 /S 

where H eq2 is given by Eq. (B.24), and r dec can be computed making use of Eq. (2.7). 

Let us now work out the perturbed energy densities. As far as matter is concerned, 
Eq. (B.10) is still valid, hence Eqs. (B.11)-(B.12) apply. As for the radiation component, one 
has 


P%- dec = pi q 2 exp [-4 (iV^-dec - N e q 2 )] (B.27) 

— P° q 2 exp [ 4 (A^,__dec AI end )] exp ( 4 (/V en( j A>.q ] ) ] exp ( 4 (iV e q ] /V,j_ mass )j 

exp [-4 (/Vo— mass - A^-dec)] exp [-4 (iV CT _ dec - N e q2 )] (B.28) 

= 3<rJ ' (g|) ' exp ^ - 2^ exp [-4 (iV 0 _ dec - N end )] . (B.29) 

By varying this expression under the constraint that A^_d ec . is fixed, one obtains 

_ 8 d In id 8 <9 In H eqi q~ eqi d In (7 eqi dN en d qn/ 

“ ~~ 3 d '0 a * 3 d(j> a * Ml 1 J 

AI 8 d 2 In H end 8 <9 2 ln H cqi ^ qi / <9 2 lnq- eqi | ^ d In cr eqi d In cr eqi \ ^ ^ d 2 N end 

0/3 3 d(j) a *d(j)/3* 3 d(f) a *d<j)p* Ml \dcj) 

(B.31) 

In these equations, the perturbations of H e q can be obtained by taking the derivatives of 
the equality between Eqs. (B.18) and (B.20). One obtains 
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At last, the perturbations of <7 eqi can be obtained from Eq. (B.19). One has 
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Let us now calculate the reheating parameter. Reheating is made of a phase of matter 
(driven by </>), inflation, matter (driven by a), radiation, and matter (driven by <f> again). 
Making use of the formulas derived above, one has 
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(B.37) 
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(B.39) 


6 V^cq,2, 

The parameter lni? rad is given by the sum of all these contributions. 

Finally, let us derive under which conditions case 2 takes place. In the limit T^ <C T CT <C 
m a < H end , one has o- eqi cr end , H 2 qi ~ m 2 a a 2 nd /(GA'I 2 ,) and pf nd ~ 3M 2 ,H 2 nd . This gives 
rise to 


H cq ~ -T CT ^exp ( 3 - 
q2 9 Mp, P V 2 M 2 , 


The condition T^ < H eq2 < To- then gives rise to 
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where W-\ is the —1 branch of the Lambert function. The lower bound is a numerical 
constant, ~ 1.69, which is close to the upper bound (B.17) \/2e 2//9 derived in case 1. The 
upper bound can be approximated in the limit T^ T^ and one obtains 


B.1.3 Case 3 



(B.42) 


In case 3, the curvaton dominates the energy content of the Universe while it is still slow 
rolling and while the inflaton has not decayed yet, but does not let the inflaton dominate 
again at later times. In order to determine the time t e q when this occurs, one needs to solve 
the equation pf q = p° q , where 
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^ ^ («end 

and, combining Eqs. (B.l) and (B.3), 
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(B.43) 


(B.44) 


The equation pt q = p° q then needs to be solved numerically. Then follows a phase of inflation 
driven by a, between cr eq and ov-mass — \/2M P1 . The realized number of e-folds is N a - mass — 
N eq = <j 2 q / (4M|j) — 1/2. During this epoch, the energy density of the inflaton field decreases 
according to 

pi- mass = Peqexp [-3 (N a - mass ~ ATeq)] = pt nd eX P ' ( R45 ) 

Then, the Universe is matter dominated until the curvaton decays. One then has 

pt *O = pL™. ( 2 = ipi- (jpY 2 ■ (B.46) 

\flg— mass/ \^g J 


From this moment, the Universe undergoes a radiation dominated era, and one then has 

/IT \ 3/2 /p \ 3/2 

= P*-*c =pL j<c (j/J . (B.47) 

One the other hand, one simply has 


P<j >-dec — 3M p [T^ P0_dec > (B.48) 

and rdec is then given by Eq. (2.7). 

Let us now work out the perturbed energy densities. As far as matter is concerned, 
Eq. (B.10) is still valid, hence Eqs. (B.11)-(B.12) apply. As for the radiation component, one 
has 


P™— dec = Pg -dec ex P [~ 4 (^fli-dec “ -^a-dec)] (B.49) 

= Pa -dec ex P I” 4 ( N <f >-dec “ N end )] exp [-4 (AT end - AT eq )] 
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By varying this expression under the constraint that N^-^ec is fixed, one obtains 
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(B.52) 


In these expressions, the derivatives of ii eq can be obtained by Taylor expanding Eqs. (B.43) 
and (B.44), which exactly match Eqs. (B.18) and (B.20). As a consequence, one obtains the 
same relations as Eqs. (B.32)-(B.33). In the same manner, the variations of cr eq follow from 
Eq. (B.44) and are given by Eqs. (B.34)-(B.35). 

Let us now calculate the reheating parameter. Reheating is made of a matter phase 
driven by 0, a phase of inflation driven by a, and a matter phase driven by a also. One has 
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The parameter lni? ra d is given by the sum of all these contributions. 

Finally, let us see under which condition case 3 takes place. In the limit T^, < T ff < 
m a < ii en d, one has cr eq ~ o- en d, #f q - 2m^^ nd /(6M^ 1 ) and pf nd ~ 3M^H^ nd . This gives 
rise to 

4-d.c = 2 (£) 7 r >. 2 »o ex p (| - riff) • < B - 56 > 

Requiring that p$_ dec < p r ^_ dec yields 


^end 
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(B.57) 


which exactly corresponds to the upper bound of case 2. 


B.2 T ff < T^ < 771(7 ^ -Rend 

In this case, the curvaton becomes massive before the inflaton decay, and it decays only 
afterwards. This corresponds to cases 4, 5 and 6. Since the last field to decay is the curvaton, 
one has cp m = a and ip T = (f>. 
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B.2.1 Case 4 

In this case, the curvaton never dominates the energy budget of the Universe. After inflation 
stops, the Universe experiences a phase of matter dominated era, followed by a phase of 
radiation. In this sense, the situation is quite similar to case 1 and these two cases share 
common formulas. Because the Universe is still dominated by the inflaton between H em j and 
r0, the oscillations of the inflaton at the bottom of its quadratic potential give rise an era of 
matter domination, and H(t) is given by Eq. (B.l). The decay of the inflaton occurs when 
H = r«£, that is at time given by 3/2(t0_d ec — tend) = 1/T0 — 1/H cnd . At this decay time, the 
energy density of the inflaton, which is instantaneously converted to radiation, is given by 

4-l.c = 1 • (B.5S) 

Then, the universe is radiation dominated and the decay of a occurs when H = T a . At this 
point, p r a _ dec is simply given by 


A-i* = pLtS- • (B.59) 

-^end 

As far as a is concerned, at the moment when it becomes massive, <7 CT _ mass is given by 
Eq. (B.4), and one has p^-mass = 3m^er^_ mass /4. Then, the curvaton behaves as a matter 
fluid in a matter dominated Universe, and one has 
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Then, the curvaton behaves as a matter fluid in a radiation dominated Universe, and one has 
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From here, rd ec is given by Eq. (2.7). 

Let us now work out the perturbed energy densities. They must be evaluated at the 
fixed background time N a _ dec , so this quantity must remain unperturbed. For the radiation 
component, one has 
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As for matter, one has 


Pa —dec — Pa —mass ex P [ 3 (N a — dec —mass)] 
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This gives rise to 
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In case 4, the reheating parameter R rad is given by the same expression as in case 1, 
namely Eq. (B.16). 

Finally, let us derive under which conditions case 4 takes place. One simply needs to 
make sure that the Universe is dominated by radiation at the decay time of the curvaton. 
Since p r a _ dec — 3M^T 2 , one has 
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B.2.2 Case 5 
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In this case, the curvaton comes to dominate again the energy content of the Universe, at 
some point between the decay of the inflaton and its own decay. First, let us notice that 
until the decay time of the inflaton, the situation is similar to case 4, so that p^_ dcc is given 
by Eq. (B.58) and p^_ dec is given by Eq. (B.61). Then, a behaves as a matter fluid in a 
radiation dominated universe, and its energy density at the time of equality between cp and 
a is given by 
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In the same manner, one has 
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Then, the Universe is dominated by a matter fluid, and one can write 
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in agreement with the fact that the curvaton dominates the energy budget at this time. On 
the other hand, the energy density of radiation is given by 
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(B.78) 


and rdec can be computed through Eq. (2.7). 

Let us now work out the perturbed energy densities. For the radiation component, 
Eq. (B.65) applies, and therefore Eqs. (B.66)-(B.67) are valid. In the same manner, for the 
matter component, Eq. (B.68) still applies, and therefore Eqs. (B.69)-(B.70) can be used. 

Let us now calculate the reheating parameter. In this case reheating is made of a matter 
phase driven by < f>, a radiation phase, and a matter phase driven by a. One has 

(R79) 

(R80) 

The parameter lni? ra d is given by the sum of these two contributions. 

Finally, let us see under which conditions on the parameters case 5 takes place. Case 
5 corresponds to the situation when equality between inflaton and curvaton energy densities 
happens between the inflaton decay and the curvaton decay, that is To- < H eq < Together 
with Eq. (B.76), this leads to 

^ 2/9 V2e 2 / 9 . (B.81) 

\ T^ / Mpi 

B.2.3 Case 6 

In this case, the curvaton dominates the energy content of the Universe from some point on 
between the moment when the inflaton becomes massive and the moment when the curvaton 
becomes massive. In order to determine the value of the Hubble parameter H cq when this 
happens, one can proceed as in case 3 and solve the equality between Eqs. (B.43) and (B.44). 
Then begins a new phase of inflation driven by a. It ends when the curvaton becomes massive, 
at Ha-m ass = m CT /\/2. The energy density contained in </>, which behaves as matter during 
this period, pf _ mass , is then given as in case 3 by Eq. (B.45). Then, the Universe is matter 
dominated and the inflaton continues to behave as matter, until it decays. One then has 

4- d« = pi-**. ( 2 = 2 pi-**. (Hj 2 . (B.82) 
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After the inflaton decay, its decay products scale as radiation and one has 
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One the other hand, one simply has 
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and r^ec can be computed through Eq. (2.7). 

Let us now work out the perturbed energy densities. For radiation, one has 
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This gives rise to 
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In these expressions, the derivatives of H eq are given by Eqs. (B.32)-(B.33) and the ones of 
<7 eq are given by Eqs. (B.34)-(B.35). As for matter, one has 
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This gives rise to 
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In case 6, the reheating parameter -R ra d is given by the same expression as in case 3, 
and Eqs. (B.53), (B.54) and (B.55) apply. 

Finally, let us see under which conditions case 6 takes place. In the limit where the cur- 
vaton remains frozen at (7 en d until equality, one has (7 eq ~ (7 en d and H eq ~ m a a en d/(V6M Pl ). 
Case 6 corresponds to the situation where this happens between the end of inflation and the 
moment when a becomes massive, i.e. m ( j/\/2 < H eq < H c n d- This leads to the condition 

\/3 « « y/6^ . (B.92) 

M P1 m a 

The lower bound y/3, is numerically very close to the upper bound derived for case 5, \/2e 2 / 9 . 


B.3 T^r <C ITlfj <C T^ < H enc [ 

In this case, the inflaton decays before the curvaton becomes massive. This corresponds to 
cases 7, 8, 9 and 10. Since the last held to decay is the curvaton, one has ip m = a and ip v = (j>. 


B.3.1 Case 7 

In this case, the curvaton is always subdominant. Between the end of inflation and the 
inflaton decay, the Universe is matter dominated and Eqs. (B.l) and (B.3) apply. Evaluated 
at the moment when cj> decays, when H = T^, this gives rise to 
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Then, the Universe is radiation dominated, and one has 
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In the limit where a is slowly rolling until the moment when it becomes massive, one then 
has 
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and Pct- mass = 3m 2 <7 2 _ mass /4. Then, a becomes massive, and one has 
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(B.95) 
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(B.98) 


One the other hand, one simply has p r a _d ec = 3M P1 T 2 — p“_d ec an d r dec is given by Eq. (2.7). 
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Let us now work out the perturbed energy densities. They must be evaluated at the 
fixed background time N a - dec , so remember that this quantity must remain unperturbed. 
For the radiation component, one has 


Pa- dec = 4-dec eX P ( N cr—dec ~ ^-dec)] 

= 4-dec eX P I -4 ( N <t ~dec “ N end)} exp [-4 (N elld - i\^_ dec )] 

= 3M^rJ ' exp [-4 (AT CT _ dec - N end )} . 

This gives rise to 
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(B.99) 
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As for the matter energy density fluctuations, one has 
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(B.104) 


This gives rise to 
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(B.106) 


In case 7, the reheating parameter R rSid is given by the same expression as in cases 1 
and 4, namely Eq. (B.16). 

Finally, let us check under which conditions the assumptions of case 7 are valid. Eq. (B.98) 
gives rise to 


Pa —dec 1 I m a ( ^end \ 


(B.107) 


Therefore, case 7 corresponds to the regime where 
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- 30 - 

























B.3.2 Case 8 


In this case the curvaton comes to dominate the energy budget of the Universe, after it be¬ 
comes massive and when the inflaton has already decayed. Until the time when the curvaton 
becomes massive, the same analysis as in case 7 apply, and Eq. (B.96) giving a n _ mass can be 
used. At the time of equality, the energy density of a is then given by 
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On the other hand, one simply has 

-4 


Pe q P ip—dec 


a, 


eq 


^(p —dec 


P 4 >—dec 


H t 


eq 




—dec 


<p 

— Pend 


H, 


eq 


H, 


end 


and solving p r eq = p° q gives rise to 
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Then, the Universe is radiation dominated and one has 
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(B.109) 

(B.110) 
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in agreement with the fact that the curvaton dominates the energy budget at this time. On 
the other hand, the energy density of radiation is given by 


r _ r 

Pa —dec Pe q 



" 2/3 , (B.113) 


and r f je C can be calculated thanks to Eq. (2.7). 

In case 8, the reheating parameter i7 ra d is given by the same expression as in case 5, 
and Eqs. (B.79) and (B.80) apply. 

Let us now work out the perturbed energy densities. For the radiation component, 
Eq. (B.101) still applies, and therefore Eqs. (B.102)-(B.103) are valid in case 8 too. For the 
matter component, Eq. (B.104) still applies, and therefore Eqs. (B.105)-(B.106) are valid in 
case 8 too. 

Finally, let us see under which conditions on the parameters case 8 takes place. Case 
8 corresponds to the situation when equality between inflaton and curvaton energy densities 
happens between the moment when the curvaton becomes massive, and the decay of the 
curvaton, that is r CT < H cq < m a /V2- Together with Eq. (B.lll), this leads to the condition 

2 5/8 e l/3 ( M V4 <<; << V2el/3 . (B.114) 

\m a J M P1 
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B.3.3 Case 9 


In this case the curvaton comes to dominate the energy budget of the Universe before it 
becomes massive but after the inflaton decay. This means that there is an extra phase of 
inflation. Until the moment when the curvaton dominates the energy content, the situation 
is similar to case 7 and Eqs. (B.93), (B.94) and (B.95) apply. This gives rise to 
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One the other hand, one has 

r _ 4> 

Pe q — P dec 


eq 


^(p —dec 


-4 


P (p—dec 


H ( 


eq 




—dec 


_ 4> 

— Pend 


H, 


eq 


H, 


end 


(B.115) 


(B.116) 


(B.117) 


The value of H eq is obtained by numerically solving p a cq = p r cq . The new phase of inflation 
terminates when ov-mass = \f2M PX and H a _ mass = m a /y/2. The realized number of e-folds is 
given by N a - mass - N eq = (o^ q - cr^_ mass )/(4M| 1 ), and one has 


Pa -mass = Peq ^ ^ ss ^j = p T eq exp [-4 {N a _ mass - 7V eq )] = p f nd exp ^2 - . 

(B.118) 

After the curvaton becomes massive, the Universe is dominated by a matter fluid until the 
curvaton decays. One then has 


r _ r 

Pa —dec Pa —mass 


P end 



(B.119) 


while p“_ dec = - p T a _ dec , and r dec is given by Eq. (2.7). 

Let us now derive the perturbed energy densities. For radiation, Eq. (B.101) is still 
valid, and Eqs. (B.102)-(B.103) still apply. As for matter, one has 
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This gives rise to 
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As explained above, H eq has to be determined numerically. However, Taylor expanding 
Eqs. (B.117) and (B.116), one can express its derivatives as 
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On the other hand, derivatives of cr eq are derived from Eq. (B.115), and read 
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Let us now calculate the reheating parameter. In case 9, reheating is made of a matter 
phase driven by 0, a radiation phase, and inflation phase driven by a and a matter phase 
driven by a. One has 


ln Cd =g ln 


In ifj-ad = - - 

rad 2 4Mp 
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(B.129) 
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The parameter In f? rad is given by the sum of these two contributions. 

Finally, let us check under which conditions case 9 takes place. In case 9, equality 
between curvaton and radiation occurs while the curvaton is still slowly rolling, but after the 
inflaton decay. In the limit T CT <C m a <CT^< ff end , one then has 


Pe q — Pend ~ <* <^end 


(B.132) 
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p -“ = = <B ' m) 

Equaling the two above relations gives rise to H cq ~ m 0 -cr en( j/(v / 6Mp 1 ). Case 9 corresponds 
to the situation where this happens between the decay of the inflaton and the moment when 
the curvaton becomes massive, that is m a /y/2 < H cq < I\. This leads to the condition 


V / 3<^<\/6^. (B.134) 

Mpi m a 

B.3.4 Case 10 

In this case equality between curvaton and inflaton energy densities occur while the inflaton 
is still massive, and when the curvaton is still slowly rolling. In order to determine the value 
of the Hubble parameter H cq when this happens, one can proceeds as in case 3 and solve the 
equality between Eqs. (B.43) and (B.44). From this moment on, the Universe experiences 
a new phase of inflation. At some point during this inflationary second epoch, the inflaton 
decays. Since, in the slow roll approximation, one has dec — m (T <70_dec/(\/6Af P1 ), this 
occurs when cr^_ dec = The number of e-folds elapsed at this point is then 

given by N^- de c — A r eq = oy q /(4M/,) — 3T^/(2m^.). Therefore, the energy density contained 
in the inflaton field when the inflaton decays is given by 
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After the inflaton decays, the Universe continues to inflate until the curvaton becomes mas¬ 
sive and Ho—mass = m a /y/2. The number of e-folds elapsed during this period is given by 
N(j -mass — iV^_dec = 3T^/(2 m^) — 1/2, so that the energy density of radiation becomes 

Pa -mass = P %-dec eX P [~ 4 (^-mass - A^-dec)] = pf* [j^j eX P ( 2 “ ■ 

(B.136) 

After the curvaton becomes massive, the Universe is matter dominated, so that the energy 
density of radiation scales as 
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(B.137) 


On the other hand, the energy density of the curvaton at this time is simply given by 
P™_dec = SM^T^ — /7_ dec , and r dec can be calculated from Eq. (2.7). 

Let us now work out the perturbed energy densities. They must be evaluated at the 
fixed background time IV 0 -_ dec , so remember that this quantity must remain unperturbed. 
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For the radiation component, one has 

Per —dec = P^-dec eX P [ — 4 (iV CT _dec — -^0-dec)] (B.138) 

= P \-dec ex P [~ 4 ( N a-dec ~ -/Vend)] exp [-4 (IV end - N eq )] 

x exp [-4 (AT eq - iV^_dec)] (B.139) 

/it \ 2/3 / 2 3 p 2 \ 

- (-£) (st - 4) exp h4 w - d “ - " md)1 ■ <B - 140) 

One can see that this expression only differs from Eq. (B.86) by an overall constant factor 
that disappears in the logarithmic derivative, and therefore, Eqs. (B.87)-(B.88) can be used. 
As for matter, Eq. (B.89) applies, and Eqs. (B.90)-(B.91) can be used. 

In case 10, the reheating parameter .R rad is given by the same expression as in cases 3 
and 6, and Eqs. (B.53), (B.54) and (B.55) apply. 

Finally, let us check under which conditions case 10 takes place. In the limit where the 
curvaton remains frozen at cr e nd until equality, one has H eq ~ m CT <7 en d/('\/6Mpi). Case 10 
corresponds to the situation where equality happens between the end of inflation and the 
decay of the inflaton, that is < H eq < i/ end . This leads to the condition 

. (B.141) 

m a Mpi m a 

C Predictions for a few models 

In this appendix, we plot the predictions of three single-field inflationary potentials when 
an extra massive scalar held a is added: Large Field inflation for which V = M 4 (</>/Af P1 ) p , 
Higgs inflation - the Starobinsky model - for which V = M 4 [l — exp( — ^2/3<f>/M Pl )] 2 , and 
Natural Inflation for which V = M 4 [ 1 + cos(</>//)]• Large Field inflation is a prototypical 
example of single-held model that predicts a value for r that is too large, Natural Inflation is 
a prototypical example that predicts a value for n s that is too small, and Higgs inflation is a 
prototypical example that predicts the right value for n s and a value for r that is sufficiently 
small. 

The value of M is set so that Eq. (2.13) hts the measured amplitude of the scalar 
power spectrum. For each potential, we investigate the 10 reheating scenarios and calculate 
n s , r and /nl- The conventions used in the following series of hgures are the following. 
In the upper left corner of each plot, the name of the model is written, where “MC” (for 
“Massive Curvaton”), a number (referring to the reheating scenario under consideration) and 
the acronym of the single-field version of the model are appended. For example, “MC7NI” 
corresponds to the natural inflation potential, where a massive scalar field is added, in the 
7 th reheating scenario. The predictions of the model are displayed as colored circles, where a 
color bar on the right indicates which parameter the color encodes. 5 For comparison purpose, 
black squares stand for the predictions of the single-field version of the model, when u> re h = 0. 
For a generic reheating scenario (that is to say, when — 1 < u) re h < 1/3), those predictions 
would span a much larger region in the observable space. Finally, the black solid lines stand 
for the marginalized joint 68% and 95% confidence level regions of the Planck 2015 data [14] 

5 Most of the time, the color encodes a end /M P i. However, when a end /M P i is not single-valued in the three 
planes (ns,r), (tisj/nl) and (/nl,)"), another quantity, cr^ nd r^/(M P 1 T a ), reflecting the relative contribution 
of cr to?(, is displayed. 
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combined with BICEP2 and Keck Array. When /nl is displayed, we simply show the 95% 
bounds for the parameters on the axes. 

In passing, let us notice that when going from one reheating scenario the other, the 
predictions evolve continuously with the models parameters. However, because all four pa¬ 
rameters m. a , T a and cr en d are varied simultaneously, this is not always obvious from only 
looking at how the color changes at the border between two adjacent reheating scenarios. 
Besides, in principle, our calculation is valid under the approximation that all events are well 
separated in time. It should therefore not be trusted close to the boundary between two 
reheating scenarios. In a Bayesian perspective in any case, this only represents a tiny region 
in parameter space and should play a negligible role in the final results. 
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C.l Large Field Inflation + Massive Scalar Field (MCLFI) 
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Figure 3. Reheating consistent slow-roll predictions for the large field inflation models with a massive 
curvaton field, when reheating scenario is of the first (left panels) and second (right panels) type. 
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Figure 4. Reheating consistent slow-roll predictions for the large field inflation models with a massive 
curvaton field, when reheating scenario is of the third (left panels) and fourth (right panels) type. 


38 










































MC6LFI 





'%P = 4 



1.5b R 

S 


-5 b R 


MC6LFI 



p—Qp —4 £ 


“00 02 04 05 OS Too 



1.5b R 


Figure 5. Reheating consistent slow-roll predictions for the large field inflation models with a massive 
curvaton field, when reheating scenario is of the fifth (left panels) and sixth (right panels) type. 
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Figure 6. Reheating consistent slow-roll predictions for the large field inflation models with a massive 
curvaton field, when reheating scenario is of the seventh (left panels) and eighth (right panels) type. 
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Figure 7. Reheating consistent slow-roll predictions for the large field inflation models with a massive 
curvaton field, when reheating scenario is of the ninth (left panels) and tenth (right panels) type. 
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C.2 Higgs Inflation + Massive Scalar Field (MCHI) 
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Figure 8. Reheating consistent slow-roll predictions for the Higgs inflation models with a massive 
curvaton held, when reheating scenario is of the first (left panels) and second (right panels) type. 
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Figure 9. Reheating consistent slow-roll predictions for the Higgs inflation models with a massive 
curvaton held, when reheating scenario is of the third (left panels) and fourth (right panels) type. 
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Figure 10. Reheating consistent slow-roll predictions for the Higgs inflation models with a massive 
curvaton field, when reheating scenario is of the fifth (left panels) and sixth (right panels) type. 
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Figure 11. Reheating consistent slow-roll predictions for the Higgs inflation models with a massive 
curvaton field, when reheating scenario is of the seventh (left panels) and eighth (right panels) type. 
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Figure 12. Reheating consistent slow-roll predictions for the Higgs inflation models with a massive 
curvaton field, when reheating scenario is of the ninth (left panels) and tenth (right panels) type. 
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C.3 Natural Inflation + Massive Scalar Field (MCNI) 
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Figure 13. Reheating consistent slow-roll predictions for the natural inflation models with a massive 
curvaton held, when reheating scenario is of the first (left panels) and second (right panels) type. 


- 47 - 









































f=HM 


0.885 0.900 0.915 0.930 0.945 0.960 0.975 0.990 

n 



1.5b" R 

d 



1.5b l 1 ^ 

d 


MC4NI 


" '• ‘"L . 

..:■* *’ 

• * * 

.. •••• 

• * * ^ .**''• *'I«S***»”* , ’ % '*** 

••• • ■* V >••'*••• 
... • ! I.v* •* ****** 


.’.’ . •••••• 

.. ... 

. ..*..••• . .%*••• 

*.**«•*.••* .•••••"* **• ••• 

• •.•••*. •••*•••••:*. •':# . 

* ■••**"** '. /«*■'! ••• 


f—3M 


f= 3M 


0.S85 0.300 07915 ff?30 0355 0950 0375 0350 005 


t-f L b 


MC4NI 


i.5b" R ^ 
d 


/=6M **’’•’•* 


10 3 f — 3M 






Ali i | 


io^ iff 1 iff 5 iff iff iff iff iff iff iff iff 


/ 

J Ml 


Figure 14. Reheating consistent slow-roll predictions for the natural inflation models with a massive 
curvaton field, when reheating scenario is of the third (left panels) and fourth (right panels) type. 
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Figure 15. Reheating consistent slow-roll predictions for the natural inflation models with a massive 
curvaton held, when reheating scenario is of the fifth (left panels) and sixth (right panels) type. 
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Figure 16. Reheating consistent slow-roll predictions for the natural inflation models with a massive 
curvaton field, when reheating scenario is of the seventh (left panels) and eighth (right panels) type. 
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Figure 17. Reheating consistent slow-roll predictions for the natural inflation models with a massive 
curvaton field, when reheating scenario is of the ninth (left panels) and tenth (right panels) type. 
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